Open environments for quantum open systems 
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The majority of quantum open system models in the literature are simplistic in the sense that 
they only explicitly account for that part of the environment that directly interacts with the system 
of interest. A quantum open system with an open environment is examined using the projection 
operator method in the weak coupling limit. The openness of environment is modelled by nonunitary 
evolution of the Lindblad form. Under certain conditions, the resulting master equation for the 
system is insensitive to the initial state of the environment and to initial entanglements between 
the system and environment for time scales greater than the environment relaxation timescales. For 
the particular case of an environment consisting of a harmonic oscillator bath, the resulting master 
equations are demonstrated to have the algebraic form for completely positive evolution. The open 
environment model is illustrated for the particular case of a system linearly coupled to an oscillator 
bath. 



I. INTRODUCTION 

There has been much interest in quantum open systems 
(particularly with applications to quantum decoherence 
[1, 2, 3] over the past several decades. [4] 

As discussed by Alicki,[4] models that have been exten- 
sively studied in the literature have arisen from one of two 
approaches. The first starts with an a priori proposal of 
a mathematical form for the quantum dynamical semi- 
group (Alicki's axiomatic approach), and proceeds with 
a study of the resulting properties of the dynamics. The 
second is to study the reduced dynamics of the system of 
interest which is part of a composite system which also 
includes environmental degrees of freedom (Alicki's con- 
structive approach). The models primarily explored for 
the constructive approach have a set of common features. 
The system of interest evolves unitarily in isolation. The 
environment (often an oscillator bath) evolves unitarily 
in isolation. The interaction between system and envi- 
ronment consists of a term in the composite hamiltonian, 
so that the composite system evolves unitarily. The ini- 
tial state of the composite system is an important feature 
of the particular model being studied and is often but not 
exclusively taken to be a factor state between an initial 
state of the system of interest and a thermal equilibrium 
state for the environment. 

The study of open systems is an attempt to describe 
how "the rest of the universe" influences the system of 
interest. However, most constructive models only include 
that part of the rest of the universe with which the rest 
of the universe directly interacts (through an interaction 
term in the composite hamiltonian) . In order to account 
for that part of the universe with which the system of in- 
terest does not directly interact, we explore properties of 
models for which the environment degrees of freedom are 
themselves open systems. In the Section II, we motivate 



the quantum models by examining a classical system and 
environment where the environment degrees of freedom 
are damped oscillators relaxing into thermal equilibrium. 
In Section III we review some important properties of the 
model we will use to simulate the open nature of the envi- 
ronment. In Section IV we proceed to use the projection 
operator method in the weak coupling regime on a com- 
posite system where the environment relaxes to thermal 
equilibrium via a nonunitary evolution of the Lindblad 
form. [5] We illustrate the results for the particular case 
of linear coupling to an oscillator bath in Section V. We 
finish in Section VI with a discussion of our results. 



II. CLASSICAL OPEN SYSTEM 

In this section we examine a classical open system 
with a damped environment to gain insight and antic- 
ipate properties of analogous quantum models. We are 
essentially following Zwanzig[6] with the addition of dis- 
sipation and thermalization to the dynamics of the envi- 
ronment. 

The system is described in terms of its coordinate x. 
The environment consists of a bath of independent oscil- 
lators with coordinates q M . The composite system equa- 
tions of motion are given by 



mx 



dU{x) 
dx 



, da 



(2.1) 



for the system coordinate, and 



+m f iuf l a fi (x) + Fp(t) 



(2.2) 
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for each environment degree of freedom. These equa- 
tions of motion correspond to the addition of dampen- 
ing (with friction coefficients -q^ = 2m fJ- "/^ ) and thermal 
noise terms to the environment portion of the closed 
system dynamics obtained from the composite system 
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Lagrangian: 

L= \mx 2 -U{x) + Y,^{Ql-^-a^x)f} (2.3) 



The Kronecker delta indicates the oscillators in the bath 
are driven by independent noise terms. The Dirac delta 
function indicates that the driving forces are white noise. 



Each oscillator is taken to be independent but driven to 
the same temperature by dissipation and white noise so 
that the following fluctuation-dissipation relation holds: 



(F„{t)F v {s)) = 4^ mfi k B T6^S(t - s). 



(2.4) 



With il M = — 7^ taken to be real (the bath 

oscillators are underdamped), the formal solutions for 
Eq. (2.2) can be written in terms of the oscillator's ini- 
tial position q^o and velocity q^o as 



1 /"' 1 

%{t)= [<^ocosfy + — (q^o +7^) sinf V] e ~ 7 '' t + / —sinn^t- s)e-^ t - s > F^ds 

a^{x(t)) - ■^-a ll {x{Q))[il ll cosQ. ll t + 7^ sin f^e -7 "' 

- /* [fi„ cos %s + 7ai sin n flS ]e-^ s da ^ t ~ s)) x(t - s)ds. 



(2.5) 



Upon substitution of Eq. (2.5) into Eq. (2.1) and rear- 
ranging, the equation of motion for the system can be 
written 



TTLX 



dU(x) f* 



dx J0 



r)(x(t),x(s);t, s)x(s)ds 
+F E1 (t) + F E2 (t) + F E3 (t) (2.6) 



where the noise terms are given by 



F E i(t) = Vm /i o;2[^ _ a/i ( a; (o))][cosO iU i+-^sinO^]e- 7 ''*, 



(2.7a) 



F E2 (t) = ^m^^ofcos^t 



sin fl^t. 



(2.7b) 



FE S (t) = f sinn^t-s)e-^ t -^F^ S )ds 



dx 



(2.7c) 



r 



and the dissipation kernel is given by 

r](x(t),x(s);t, s) = ^m M u^{[cosf^(t - s) 

+^ (t _ a)]e -^t-s)da^x(s))da^(t)) 

7^ ox ox 

In the past, the properties of noise terms like F E1 {t) and 



F E 2 (t) have generally been extracted from the statistical 
distributions of the initial conditions. With the introduc- 
tion of dissipation in the environment, those noise terms 
can readily be seen to be transient terms on the time 
scales of the environment (as determined by 7^). Thus 
the details of the initial state of the environment are not 
important to the long term system dynamics. 
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The correlations (FE3(t)FE3(T)} of the remaining noise 
term depend upon the correlations of the individual os- 
cillators' noise terms via Eq. (2.4), and can be written 

(F E3 {x,t)F E3 {x',r)) = k B Tr)(x,x';t,T) 
+j j hl cosn n( t + T ) - u l cosQ^-r) 

^^iv^^^ yy . (2.9) 

The first term on the right hand side of Eq. (2.9) is the 
long term correlation function of the effective noise. The 
second term contains a factor of e~ 7 ^* +T \ and thus is 
transient for long timescales. 

Since t]{x, x'; t, r) is a narrow function of t — t, we can 
make a Markov approximation in Eq. (2.6): 

dU(x) 



mx = — - 



dx 



- tj(x,x)x + F s (x,t), 



where 



r](x,x') = ^2to 



da^x)) da^x') 



dx 



dx 1 



(2.10) 



(2.11) 



The corresponding fluctuation-dissipation relation is 

(F s (x,t)F s (x',T)) = 2k B T V (x,x')8(t-T). (2.12) 

The role of the spatial correlations of the noise in quan- 
tum decoherence and the lack of importance of those 
correlations in classical phenomena has been discussed 
elsewhere. [7, 8] 

Thus, adding "fast" thermal relaxation to the environ- 
ment of the system of interest lead to Markovian equa- 
tions of motion with the usual fluctuation-dissipation re- 
lations. Memory effects due to the details of the initial 
environment state (including possible correlations with 
the initial system state) are erased on the environment 
relaxation time scales. This result provides excellent mo- 
tivation for exploring similar models in the quantum me- 
chanical domain. 



III. PROPERTIES OF QUANTUM OPEN 
ENVIRONMENT 



As with the classical model, in order to account for the 
open nature of the environment we need to incorporate 
modifications to the dynamics of the environment degrees 
of freedom. The environment is modelled as a set of 
independent oscillators whose evolution is governed by a 
Markovian master equation: 



(3.1) 



The generator L of the evolution is taken to be of the 
Lindblad form [5]: 

d P t\ l 

at = L[p] 



Formally, the solution to Eq. (3.2) is 
p(t) = A(t)[p] = e Lt [p], 



(3.3) 



expectation values are defined via the trace operation, 
allowing the definition of the adjoint representation of 
the evolution operator A* (t) [O] via 



Tr[pA*(t)[0}}=Tr[A(t)[p]0], 
and for the generator L* 

Tr[pL*[0}}=Tr[L[p}0}. 



(3.4) 



(3.5) 



The adjoint representation of L (i.e. Heisenberg picture) 
is given by 

L *[0] = - 1 [H,0}+ ^ £^[0,Vj+[^,0]^. (3.6) 



To model the relaxation of the environment, we will 
use a subset of a family of master equations that have 
been studied extensively in the literature [9, 10, 11, 12] 
where the Lindblad form of the master equation can be 
rewritten as: 

I - >o,p]-4(A + ,)fe,W}] 

+^-ri\p,ti,p}]-^ L h hp]] 

D, 



+ £^([ q ,[p,p}] + [p, [ q ,p]]). 



(3.7) 



The details of the model are determined by the specifi- 
cation of the diffusion coefficients D qq , D pp and D pq and 
damping constants fi and A, subject to the constraints: [9] 



D qq > 
D pp > 



D D - D 2 > (—) 2 
±j qq ±j pp u pq i I, 2 J ■ 



The nominal hamiltonian is 



p 2 mwg 2 
2m 2 



(3.8) 



(3.9) 



p is taken to be the (unique) asymptotic state of A(t), 
that is 



Aoo[p] = Km A(t)[p] = p 

t—>OQ 



(3.10) 



for any initial environment state p. Since p is a stationary 
state of L 



L[p] = 0. 



(3.11) 



The asymptotic behavior of the environment is encapsu- 
A 00 [p]=pW ( 3.i2) 



lated in A^, so that 
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for any state p. Although Aoo is defined only on the 
space of density operators (including pure states), we will 
need to extend its domain. If {|efe)} is a basis for the 
environment Hilbcrt space, then 



Aoo[|efe)(efe|] = p 
for any k. Furthermore, if 

\i>nm) = "^f d e «) + l e ™)) 

\Xnm) = ^|(l e «) +i\ e m)) 

then 

hoo[\4>nm)(lpnm\] = P 



(3.13) 



(3.14) 



= ^(^cc[\en}(e n \} + A 00 [|e m )(e m |] 
+Aoo[|e n )(e m |] + Aoc[\e m )(e n \]) 



= P+ ^(Aoc[\e n )(e m \} 
+A 00 [|e m )(e„|]), 



(3.15) 



and 

Aqo [|Xnm) {Xr, 



2( A oo[|e„)(e„|] + A 00 [|e m )(e m |] 
+iA 00 [|e„)(e m |] - ik 00 [\e m ){e n \ 



P + «-(A 00 [|e„)(e J 



(3.16) 



(3.17) 



-iA 0O [|e m )(e„|]), 

Which implies 

Aoo[|e„)(e m |] = S nm p. 
For an arbitrary operator O 

Aoo[0] = ^O nm A 00 [|e„)(e m |] 

nm 

^ mm p 

m 

= Tr[0]p (3.18) 

There is in general no guarantee that an asymptotic 
state exists for arbitrary L, that is, not all choices for 
the parameters D qq , D ppi etc. will be appropriate to 
model the a system dynamically relaxing to equilibrium. 
For example, if the parameters satisfy 



_ A + a hjJ 
Mod = — : — ftmwcoth 



2fc B T' 



A — a h fvjj 

lJ aa — COth 



2k B T' 



(3.19a) 



(3.19b) 



D pq = 0, 



(3.19c) 



then the Gibbs state is the asymptotic state. [9, 10] On 
the other hand, if 



D n 



JDpp — 



2mfl 1 



HXmuj 2 

2n ' 



n hXfl 
" 2ft ' 



(3.20a) 



(3.20b) 



(3.20c) 



(where ft 2 = u 2 — p 2 ) there can be persistent pure 
states. [11, 12] Since we are considering the primary ef- 
fect of the openness of the environment to be effective- 
ness which relax the environment towards an equilibrium 
state, we would only consider those choices of parameters 
for which there is a unique asymptotic state. 

Sandulescu and Scutaru have determined the time de- 
pendence and asymptotic behavior of various moments 
of p and q, which will be useful in later calculations. The 
evolution of the first order moments is given by 

<M = -(A -/*)<«,) + l<p) (3.21a) 



dt 



= -(X + p)(p) -muj 2 {p), (3.21b) 



which illustrates the role of the dissipation coefficients A 
and p. With 



(3.22) 

the evolution of the second order moments is given by 

d( P 2 ) 



dt 



2mu) 2 {S pq ) - 2(A + p)(p 2 ) + 2Dpp (3.23a) 



^£t = ^(S pq )-2(X-p J )(q 2 )+2D qq (3.23b) 

- Up 2 ) - mu 2 (q 2 ) - 2X(S pq ) + 2Dpq (3.23c) 

which illustrates the role of the diffusion coefficients D pp , 
D qq and D pq . 

The asymptotic first order moments are 

(p>oo = (3.24a) 



■5 



(<Z)co = 0. 

For the second order moments we have 

1 



(3.24b) 



(P 2 )c 



2A(A 2 



or 



V) 



{m 2 LU 4 D qq 



-[2A(A - fj.) + lu 2 ]D p 
-2mw 2 [A — fj]D pq }, 



(3.25) 



{q 2 )oo = 



2(mw) 2 A(A 2 +uj 2 -m 2 ) 
+w 2 D pp + 2mw 2 (A + /u)-D pg ] 



[(m«) 2 u 2 D w 



(3.26) 



and 

{Spq} oo 



2(mA(A 2 +^- M2 ) h(A + "X™")^" 



+(A - Ai)£» PP + 2m(A 2 - M 2 )L>p 



(3.27) 



We will also need the time dependence of q(t) in the 
Heisenberg picture. This can be readily extracted from 
Sandulcscu and Scutaru 's results for the time depen- 
dence for the first order moments. [9] The result is 



q H (t) = A(t)*[q] 
= q(cosQt 



ft. 



+ 



P 



sin nte~ xt 



sin Vlt)e 



with 



n 



(3.28) 



(3.29) 



and ft taken to be real (the oscillators are underdamped). 

The specific environment we will employ consists of a 
set of independent oscillators, each subject to evolution 
of the form Eq. (3.7), with position operators {qi} and 
associated parameters {fit}, {K}, etc. Thus the Hilbert 
space of the environment is actually the tensor product 
of the Hilbert space corresponding to each environment 
degree of freedom. The interaction between the environ- 
ment and system is accounted for by adding an interac- 
tion term to the composite system's hamiltonian: 



(E) 



(3.30) 



We will need the correlation functions (q n (t)q m } in the 
next section. Using Eq. (3.28), these correlations can be 
written in terms of the asymptotic correlations as 



(«n(*)«m> = Tr[A* (t)[q n }q m p\ 



= s n 



l [(q 2 ) oc (cosm + ^sinf7t)e At 



+ 



^{Spqn) c 



ih 



2mfl 



sin fit e~ Xt ], (3.31) 



where we have made use of the relation 

PnQn = \{{Pn,q n } + [Pn,q n ]) 



J pqn 



ih 



(3.32) 



IV. 



QUANTUM MASTER EQUATION FROM 
AN OPEN ENVIRONMENT 



To construct the master equation, we will use the pro- 
jection operator method in the weak coupling regime, 
largely following the presentation of Alicki and Lendi.[13] 
The composite system (the system of interest plus its en- 
vironment) is taken to evolve according to Eq. (3.1). The 
generator of this evolution is the combination of the nom- 
inal dynamics and an interaction: 



Lj 



(4.1) 



where the system and environment dynamics each con- 
tribute separately: 



Let — Lc 



L, 



(4.2) 



The generators for the system and environment nominal 
evolutions act on on the corresponding subspaces, so that 



LslA^®B^]=L { s s} lA^}®B( E \ (4.3) 



and 



L E [A^ ® = ® L { p[B^] (4.4) 
for all A^ and B^ E \ is taken to be unitary, and 

( E} 

L E is taken to be of the Lindblad form. It is useful to 
note that necessarily L$ and L E commute. 

The projection operator is defined in terms of a partial 
trace: 



P [O]=Tv E [O]®~p {E \ 



(4.5) 



We wish to establish some important relations between 
Po and the generators Ls and L E . 

L E P [O] = L E [Ti E [0]®p^] 

= Tr E [O]®4 s) [p( E )]=0 (4.6) 

( E} 

for any operator O so that L E P = 0. L E generates 
trace preserving evolution so that 



Tr B [4 B) [O( £ )]] = 0. 
for all {E \ Thus 
P L E [A^ ® B^} 



P [A< 



= (Tr E [L^[B^]])A^ 
= 0, 



(4.7) 



(4i 
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for all and B^ so that in general P L E = 0. 

= L { s S) [A^}^p^Tr^[B^} 
= L s [[AW]®pWTrW[BW}} 
= L s Po[A^®bW} (4.9) 

for all A (S ) and B^ so that in general PqL s = L s Po- 
Defining a second projection Pi = 1 — Po, it is easy to 
see that P\Ls = LsP\. In order to focus on the re- 
duced dynamics of the system of interest, we study the 
dynamics of Pop in certain approximations. Applying 
the projectors to the master equation for the composite 
system and making use of the idempotent property of 
projections produces the following equations: 



dP p 

dt 



dPip 
dt 



= PoLP PoP + P LP 1 P 1 p, 



P 1 LP 1 P lP + PxLPoPop. 



(4.10a) 



(4.10b) 



Eq. (4.10b) can be formally integrated and substituted 
into Eq. (4.10a), so that 



OPpp 
dt 



= P LP P p{t) 

+P o LP 1 [e p ^ t P 1 p(0) 

rt 



,P 1 LP l (t-a) 



P^PoPopis^ds. (4.11) 



Using the relations between the projectors and generators 
discussed above, we have 



PqLPi — PqLiPi, 
P x LPo = PiLjPo, 
P Q LP Q = L S P Q + P Q L I P Q , 

so that Eq. (4.11) can be written 

dPoP 



(4.12) 



dt 



= (L s + P Q L I P )P p(t) + P i/-Pie M Pip(0) 

+ / P L I P 1 e L ^P 1 L I P P p{s)ds. (4.13) 
Jo 



This result is exact, within the constraints placed on the 
model so far. 

The interaction, as specified by Eq. (3.30), can be writ- 
ten 



1 



Lip] = jjj[Ui,0] 



This form, along with the first order moments of the 
environment asymptotic state, simplifies the first term 
of the right hand side of Eq. (4.13). Specifically, for any 
operator O, 



PvLjPolO} = Tr B [l J2^ S) ® <?i £) >Tr £ [0] ® p^}} ® p 



n( E ) 



■p( E \ 



(4.15) 



r 



From Eq. (3.24) each term in the sum has a factor so that Eq. (4.13) can be written 
(Qrf^)oo — so that 



PoLjPq = 0. 



Furthermore 



and 



P L I P 1 =P L I (1-P ) = P Lj 



P 1 L I P Q = (1 - P^LjPo = LiP , 



(4.16) 



(4.17) 



(4.16 



OPpp 
dt 



L s Pop(t) 
+ PoLjP^P^iO) 

+ f PoLje^-^LjPvPop^ds]. (4.19) 
Jo 

The second term on the right hand side of Eq. (4.19) 

PoLjP^P^iQ) (4.20) 

represents a transient term which depends upon the ini- 
tial state of the composite system. For most constructive 
models the initial state is taken to be a factored state 
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of an arbitrary system state and the environment in its 
asymptotic state: 

p(0)=p( s >(0)®/5^. (4.21) 

For this type of initial condition p(0) = Pop(Q) so that 

Pip(O) = 0. (4.22) 

However, because we have added relaxation to the dy- 
namics of the environment, this type of factoring as- 
sumption is not necessary. We will be taking the weak 
coupling limit and so we will show that Eq. (4.20) is (ap- 
proximately) zero to the lowest nonvanishing order of the 
interaction in the remaining terms of Eq. (4.13), which 
turns out to be second order. Formally we can write 



e Lt = e L t 



f e L{t - s) L ie LoS ds. (4.23) 
Jo 



Substituting Eq. (4.23) into Eq. (4.20) and keeping only 
second order, we have the approximation 

P L I P 1 e Lt P lP {Q) w PoLjP.e^P.piO) 

+P L I P 1 [ e L °(*-^i/e LoS rfsPip(0). (4.24) 
Jo 

We are interested in timescales t which are assumed to 
be much longer than the relaxation timescales of the en- 
vironment, so that 



,L s t A 



Eoo- 



(4.25) 



Similarly, for the integral from t to s, either t or t — s 
(or both) is large compared to the relaxation timescales 
of the environment, so that cither 



e L (t-s) e L s (t-s) A 



Eoo 



or 



e LoS w e Lss A 



Eoo 



(4.26) 



(4.27) 



or both. From Eq. (3.18), for an operator O decomposed 
in terms of the basis for the environment {|efc)} and the 
basis for the system {l^v)} we have 

Aboo[0] = O am< i3 n A Eoo [\(fi a ) <g> \e m )(<j) a \ <g> (e m \] 
= O amt p„\<)) a }(4> a \ <g> A^,[|e TO )(e TO |] 
= O am ^ n \(j) a ){4> a \ <g> <5 m „/5 (B) 



Tv E [0} 



PoO. 



(4.28) 



When Eq. (4.28) is applied to Eq. (4.24) with cither 
Eq. (4.28) or Eq. (4.28) applying, then all terms in 
Eq. (4.24) end up with factors of either PqP\ or P\Pq 1 
both of which are 0. Using this result and keeping only 



up to second order in the interaction term, the master 
equation Eq. (4.19) becomes 



dPpp 
dt 



L S Pop(t) 



+ f P L I e L ^ t -^L I P P p{s)ds. (4.29) 
Jo 

We can rewrite Eq. (4.29)using Eq. (4.17) and 
Eq. (4.18) to get 



dPpp 
dt 



L s Pop{t) 



+ 



f P L I P 1 e L ^ t -^P 1 L I P Q P Q p( S )ds(4.30) 
Jo 



In the integrand in Eq. (4.30) we see the factor 
p 1 gio(t-s)p 1 w hi c ;h (following the discussion above) will 
be zero for t — s longer than environment relaxation 
timescales so that the primary contribution to the in- 
tegral is for s t. This yields the simplest Markovian 
master equation we will extract from the model: 



dP oP 
dt 



= L s P oP (t) 

+P L I L I P oP (t). 



(4.31) 



Using the particular form of the interaction given by 
Eq. (4.14) and noting that the independence of the oscil- 
lators comprising the environment implies 

(q n qm)oo = 8 nm (q*) oo (4.32) 

allows us to rewrite Eq. (4.14) in a final form: 

dPoP 



dt 



LsPopit) 



-|r£<?X[K,[K,PoP(i)]]- (4.33) 



Eq. (4.33) represents the added effect of environment in- 
duced noise on the system's dynamics which is responsi- 
ble for phenomena such as quantum decoherence. How- 
ever, additional effects such as dissipation are not present 
and will require a more careful handling of the Markov 
approximations. 

To reconsider the Markov approximations, we return 
to Eq. (4.29). The unitarity of the isolated system's evo- 
lution implies that 



e Lst [AB] = e Lst [A]e Lst [B]. 



(4.34) 



The naive Markov approximation is introduced into 
Eq. (4.29) by examining 

e £o(t-«) Ll p op(s ) = e LE ^- s) e Ls ^L I P p{s) 



L E (t-s)L s (t-s) 



[Ui,P p(s)} 



= l e L ^(*- s )[(e i ^*- s )[/ 7 ), (e Ls ^ t ~ s '>P Q p(s))} 



L E (t-s)^ e L s (t-s) 



ih 



U!),Pop(t)] (4.35) 
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The Markovian master equation can be written as 
dPoP 



dt 



L S Pop(t) + KP p{t) 



where 



K[P p] 



PoLje^LjPoPopds, 



(4.36) 



(4.37) 



which becomes for our particular model: 

K[P p(t)] = - ¥ J P o [tfi,e L *W[(e is M0}),P o p(t)]d« 



1 r°° 

/ d S {J2Tr E [[V,^ ®qi E \e L ^l(e L ^(V^ ® q {E) )), Pop(t)}}} ® P iE) }. (4.38) 
n Jo 



r 



Since Le acts only on the environment, we can make use 
of the cyclic property of the trace over the environment 
to rewrite Eq. (4.38) as 



K[P p(t)} = 



1 f°° 

-Jf ds{J2Tr E [[vW® q i E \s), 

nm 

[Vi S) (s)0qi E \Pop(t)}}}®p^}, (4.39) 



where 



(4.40b) 



In terms of the correlation functions specified by 
Eq. (3.31) we can write 



K[P p(t)} 



1 r°° 



(q n (s)q n )V n V n (s)Pop(t) 
- (q n qn(s))V n P p{t)V n (s) 
- (qn(s)q n }V n (s)P Q p(t)V n 
+ (q n q n (s))Pop{t)V n (s)V n }, 



(4.41) 



where V n = Vn S ^ ® 1^ and 1^ is simply the identify 
operator on the environment subspace. We note that here 
V n (s) — e Lss V n is not the Heisenberg evolved operator, 
that is V n (s) ^ e L 's s V n . 

While Eq. (4.39) is similar to previous results, the main 
difference comes from the nature of L^\t). Previous 

authors have taken L E (t) to necessarily generate uni- 
tary evolution in order to use properties analogous to 
Eq. (4.34) which does not generally apply to nonunitary 
evolution. 



Since Eq. (4.39) has the same mathematical structure 
as results obtained without dissipative effects in the envi- 
ronment, we expect similar shortcomings. In particular, 
to insure complete positivity for the reduced dynamics, 
we will now apply an averaging process (sometimes ref- 
ereed to as the Rotating Wave Approximation), [13, 14] 
defined by: 



K = lim - 

a— too a 



e- LoT Ke LoT dT. 



(4.42) 



(4.40a) Since L E P 



(4.43) 



so that with Eq. (4.2) and Eq. (4.37) we can write 
1 



K = lim — I e 

a^oo a J 



LsT Ke LsT dr. (AAA) 
Using Eq. (4.34) and Eq. (4.38), Eq. (4.44) becomes 

K[P p(t)} = - ¥ }™-J o drl ds^{ 

(qn(s)q n )V n (-T)V n (s - r)P oP {t) 

- (q n qn(s))V n {-T)P p{t)V n (s - r) 

- (q n {s)q n )V n {s - T)P p(t)V n (-T) 

+ (q n qn(s))PoP(t)V n (s - r)K(-r)|4.45) 

The operators {V n } can be decomposed in terms of the 
energy eigenstates of the system hamiltonian -ff(s): 



E E 

E^>- 



\p)(p\V n \v)(v\ 



Aoj e u — e l ,=HAijj 
.Aw 



(4.46) 



Aw 
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It is clear that with this decomposition 

Vn,-A» = Vl Au . (4.47) 
The time dependence of the operators {V^(s)} becomes 
V n (s) = J2e- iA » s V nAu . (4.48) 

When Eq. (4.48) is substituted into Eq. (4.45), there will 
be oscillating terms which the integral over r will cancel, 
via 

lim - / dTe^' +A ^ = <W,-Au, (4.49) 

a^oo a Jq 

so that Eq. (4.45) becomes 

K[P 0P (t)} = - p / dsJ2 ^ AUJS { 
(q n (s)q n ) 

\QnQn 

-(qn(s)q n )Vl AuJ P p(t)V nAu > 

+ (qnqn(s))Pop(t)Vl AbJ V nAu] }.(A.50) 

I 



and 



where 



The denominator in the right hand side of Eq. (4.53) is 
the product of two sums of squares, and hence is guar- 
anteed to be positive. The numerator is quadratic in Au> 



With the definition 

/>oo 

h„,Au, + tS'n.Ao, = / dse- lA ^ s (q n (s)q n ), (4.51) 

Jo 

and using Eq. (4.47) we can write Eq. (4.50) as 
K[P p(t)} = ±_S nA „V nAu Vl Au ,P p(t)} 

n,Acj 

+p ]T h nA „([vl Au P p(t),V nAoj ] 

n,Auj 

+K Au ,P p(t)V nAu }). (4.52) 

The first term in the right hand side of Eq. (4.52) is 
simply an additional hamiltonian term. The remaining 
terms (those with h nAuJ ) are in the Lindblad form if 
h nA w is positive. Using Eq. (3.31) and the moments 
given in Eq. (3.26) and Eq. (3.27) in Eq. (4.51) we find 



(4.53) 



(4.54) 



(4.55a) 
(4.55b) 
(4.55c) 
(4.55d) 

I 

of the form 

y = ciAlu 2 + bAuj + c 

= {m 2 n D qqn )Auj 2 - m n X n H Au 
+m 2 n (X n + p n ) 2 D qqn + D 

ppn + 

2(A„ + p n )m n D pqn . (4.56) 



h nA uj 



[(A„ + p n ) 2 + Au} 2 ]m 2 n D qqn + Dppn + m n \ n hAoj + 2(A» + p n )m n D p 

[\l + (n n + Au) 2 ][x 2 n + (n n -Au) 2 } 



~>n A u> — 



C an + C ln Auj + C 2 „Am 2 + C 3n Am 3 
2Am 2 [A„ + (f>„ + Aw)2][A2 + (fi„ - Alu) 2 } ' 



C Qn = hm n X n (X 2 n + n n ) 2 , 



Ci„ = mi - 'SX 2 n )(p + X) 2 + (A„ + n 2 n ) 2 ]m 2 n D qqn + (ft„ - 3A„)(2( Mn + X n )m n D pqn + D ppn ) 



C 2n = -hXm n {X 2 n +tt 2 n ), 



C 3n = -m n D qqn [X n + fl n + (p n + X n ) } - 2 D pqn m n {p n + X n ) - D p 
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If y is positive, then h n ^ u is positive as well. The coeffi- 
cient of the quadratic term is positive, thus y has positive 
concavity. If y = has no real roots, then y must be pos- 
itive, which can be tested by the condition 4ac — b 2 > 0. 
Upon some rearrangement, we can write 



Aac - b 2 



Am n {[m n D qqn (X 



+D qqn D ppn -D 2 pqn -(^) 2 }. (4.57) 



The first term inside the braces is a square and hence 
positive, while the remaining terms satisfy Eq. (3.8) and 
so the resulting expression is necessarily positive. Thus 
it is sufficient that our open system model for the oscilla- 
tors is of the Lindblad form (as discussed in the previous 
section) to guarantee that the rotating wave approxima- 
tion of the weak coupling limit generates a Lindblad form 
for the evolution of the system. 



V. EXAMPLE: OSCILLATOR LINEARLY 
COUPLED TO BATH 

In this section we illustrate our results with a test 
model consisting of an oscillator linearly coupled to a 
damped oscillator bath. To simplify notation we take p 
to be the reduced density operator for the system. The 
system has a nominal hamiltonian given by 



Hs = \m s u 2 Q 2 + JL, 



(5-1) 



and the interaction with the environment is given by 

U! = Y, C nQ®q n - (5.2) 



With this choice, we can write the operators V^Aw m 
terms of the system creation and annihilation operators 
a) and a 



V n ,Acj — C n 



{a <5a )• (5-3) 



The contribution to the hamiltonian through K, as it 
appears in Eq. (4.52), becomes 

AH = - ^ (S , „^AwK ; AwU„ ! Aw) t 



n,Atj 



Mw s (aa t - -) + AE, 



(5-4) 



where 



C 2 



) (5-5) 



and 



5E : 



y 

2m S Lu s 



(5.6) 



Thus Slos is simply a frequency shift for the system and 
SE is a C-number shift in the energy. The remaining 
contributions to K essentially are of the form of Eq. (3.7) 
with fi = 0, that is 



K[p] = - l -[M>u s {a<J -i) + AE,p] 

-^Hh{p,p}} - \p, U,p}}) 

-£EE [q , [q ,p]]-^l\p,\p,p]]. 



(5.7) 



The dissipation and diffusion coefficients are given by 



A 

Dpp 



E 



2hmsuJs 

IL 

(m s uJ S ) 2 Dqq 



{h r , 



(5.8) 



If the environment degrees of freedom are thermal- 
ized (i.e. driven asymptotically to a Gibbs state) then 
Eq. (3.19) holds and the dissipation and diffusion coeffi- 
cients are given by 
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m 2 (A 2 + (fi„ - w s ) 2 )(A 2 + (fi„ + ws ) 2 ) 
_ C 2 ficoth(^)[(A„ + /i„)(A 2 + » 2 ) + (Anj-^Vlj 

and the frequency shift is given by 

S 2m 2 Ws (A 2 + (fi„ - c s ) 2 )(A 2 + (fi„ + ^ s ) 2 ) 1 ' ' 



From Eq. (3.21), the dissipative time scales of the en- 
vironment are determined by A n — \i n and A„ + n n . If 
dissipation is weak, then 



and 



A n — fJ-ni A„ - /l n « U„ 



(5.11) 



(5.12) 



The dissipation and diffusion coefficients each have a fac- 
tor of 

1 



1 



{w n - uj s ) 2 



(5.13) 



from which we see that for a weakly damped environ- 
ment, the greatest effect is from the oscillators in the 
environment with frequencies close to the system's fre- 
quency. With this approximation, we can further sim- 
plify the dissipation and diffusion coefficients with: 



A 



E 



c 2 



4m 2 wf 7 A„ 



hm s u s fiws \ - 

Dpp w coth — — — > 



2k B T 4m|w|A n 



r 2 

5^n-(5.i4) 



Thus the parameters of K satisfy Eq. (3.19) and the sys- 
tem is driven towards thermal equilibrium by an effective 
evolution of the Lindblad form. 



VI. CONCLUSIONS 

We have constructed open environment models to ac- 
count for the environment's interaction with the "rest of 
the universe" . In Section II, our brief investigation of a 
classical model provides some insight and expectations 
for the development of a quantum mechanical model. 
The review the open system model in Section III served 
to establish many properties used in our derivation of the 
effective master equation in Section IV. Although Sec- 
tion IV largely follows previous work, the introduction of 
nonunitary evolution for the environment provided some 
novel aspects to the master equation derivation. Wc 
were able to show that an effective master equation of 
the Lindblad form could be obtained for a rich family of 
dissipative environment models. Our illustration of the 
resulting master equation with a bilinear environment- 
system interaction provided an demonstration in which 
an environment dynamically driven toward thermal equi- 
librium can naturally result in the dynamical thcrmaliza- 
tion of the system of interest. 
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